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ABSTRACT 


A non-iterative quadrature algorithm is presented which 

is based on a cubical spline interpolant. Since the underlying 

interpolant can even interpolate data which is unequally spaced 

in abscissa, the proposed quadrature algorithm can equally 

handle this type of data or equispaced data. Since the under- 
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lying interpolant is the smoothest C function which interpolates 
all the data, numerical quadrature by this technique minimizes 
spurious contributions to the integral which often result by 
least squares techniques, Simpson's 1/3 or other similar multiple 
r ules which implicitly employ a polynomial interpolate. A direct 
comparison for equispaced data between this algorithm, Simpson's 
1/3 rule, and the exact analytic quadrature for a wide range of 
polynomial and transcendental integrands shows that the algorithm 
is generally five times more accurate than Simpson's 1/3 
rule . 

The algorithm requires as input information i) the data to 
be integrated and ii) some approximation to the leading and 
trailing slopes of the interpolant. For the cases of analytical 
integrands condition ii) is easily met by the elementary calculus. 
For tabular data, the loosely coupled nature of the underlying 
interpolant is shown to manifest itself by confining inaccurate 
leading and trailing slopes effects on the interpolant to at 


worst the first and last three data points. By sampling 
outside the interval of interest one can then obtain more accurate, 
smooth interpolation and therefore quadrature on the inner 
interval of the interpolant 's domain. 


INTRODUCTION 


Determining integrals of tabulated empirical, or otherwise 
discretely known data is a general numerical problem which is 
often encountered in scientific work. This task is often 
complicated by the constraints that result from data sampling . 
Rarely is one fortunate in having the data equally spaced, or, 
if one is so fortunate, the number of data points is generally 
not some exact multiple o:f specified numbers, N, peculiar to 
the ordinary multiple rules such as those of Simpson, Hardy, 
or Weddle. Even in the case when the data is equispaced and 
such a multiple of N, the underlying interpolation of the 
discrete data is by a continuous function which is a piecewise 
polynomial of order N-l, just as one would obtain by Lagrange 
interpolation over N points. The pitfalls of Lagrange inter- 
polation are well known and are generally to be avoided in 
numerical work, especially as N becomes large. The basic problem 
with Lagrange interpolation is that it oscillates through the 
domain of interpolation in order to interpolate every data 
point. + 

If the underlying objective in taking discrete data and 
interpolating the results is to approximate a continuous function 


+ For graphic examples cf. Thompson p. 18. 


Lagrange interpolauts may introduce much spurious content. 

The integral of such a poor approximation of the continuous 
function may therefore be a poor approximation to the integral 
of the underlying continuous function. Least squaring the 
discrete approximations of the integrand to some functional 
form is also used for numerical quadratures of equally or un- 
equally spaced data. This technique is not useful if one is 
trying to infer from the integration something about the under- 
lying functional form of the data. This is often the purpose 
in taking moments of empirical distributions. 

An alternative quadrature technique which reduces the 
problems of injecting spurious information in the interpolation 
step prior to quadrature will now be discussed. Often it is 
an implicit or explicit assumption in sampling with discrete 
counters (generic) that the underlying structure of ’’events” 

evolves smoothly in time, increasing energy or other variable. + 

» 

An example (explicit) of such a "counter” would be the discrete 
monotonic tabulation of a continuously differentiable function 
at a number of points (= "events”) in its domain. Another example 
(implicit) would be a series of geiger counters with different 
energy thresholds monitoring the distribution of energetic 
solar electrons. The assumption is rather common. When discrete 
data are taken under this assumption (and all instrumental noises 

♦This assumption will be referred to as the smooth distribution 
hypothesis 


are removed) the smoothest curve which interpolates the data 
is the only interpolant which does not interject further 
assumptions about the underlying distribution of "events". 
Theoretical cubical splines are in the class of such 
interpolants . 

Recently cubical spline interpolative techniques have 

been reduced to non-iterative algorithmic forms of exceptional 

computational stability and accuracy suitable for digital 

computers (Thompson 1970, 1971). The resulting interpolant 

S(x) interpolates every data point regardless of spacing in 

such a manner that S(x) possesses a maximal smoothness property 

2 

as compared with any other C function which interpolates all 
the data, (Holladay, 1957) • By construction, the interpolant 
on any interval between data points "knows" about the variation 
of the data in the neighboring intervals so that it may make 
the "smoothest" transition from that interval to the next and 
still interpolate; all the subsequent data . This contrasts 
sharply with the local nature of the interpolants of the multiple 
rules whose derivatives are discontinuous at data points indexed 
jN, {j = 0,1,2, . 

The present quadrature algorithm exploits the piecewise 
(in general different) cubic functional form of S(x) on each 
subinterval between data points to do the necessary integrals 
via Simpson's 1/3 rule which is exact for polynomials of order 
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- 3. Since the spline interpolant considered over the entire 


interval of interest is one of maximal "smoothness " , sharp 
junctures in the interpolant are avoided. 

general scheme 


The integration technique has been formulated as a com- 
puter subroutine, INTEG, which computes numerical integrals 
for tabulated data {(x^ U (x^) , i - 1, NPTS - number of data 
points) equally or unequally spaced in abscissa. This modular 
routine is designed to be compatible and interface with and 
use the routines developed for spline interpolation by Thompson 
(1970). The spline function S(x) is fitted to the data in the 
non-iterative method of Thompson (1970) using the computer 
subroutine SPLN2.+ 

S(x) is continuous, and therefore 


A 

/ 

x 


NPTS b NPTS-1 

E 

i=l 


U(x)dx 


b 

•'a 


S(x)dx - 


A 

/ 


i+1 


S(x)dx , 


( 1 ) 


where U(x) is the underlying (smooth) distribution from which 
samples are taken. Employing the interpolative power Which 
knowledge of S(x) implies, and using the fact that S(x) restricted 


^Copies of this and other routines referenced in Thompson 
1970 may be found in Appendix B. 


i£M -*»***• ftp* 


to [x^, X i + x^ is a cubic polynomial for which Simpson's 1/3 
rule is exact, we obtain: 


S(x)dx * — (S(x i ) + 4S(x jL + 6 i > + S ^ x i+ i^ 


where 


X. _ - X, 

i+1 i 


Since S(x i > ■ U(x i > 


NPTS-1 


U(x)dx 


- | ^ 6 i (0(x i ) + 4S( Xi + 6 ± ) + U(x i+1 )) 


COMPARISON OF SIMPSON AND SPLINE QUADRATURE 

— . -t * — 

It is worthwhile to inquire whether this technique has any 
computational advantages over the equispaced rules. To 
demonstrate the advantage of spline quadrature we consider the 
following class of numerical problems: integrals without ante- 

derivatives. Although these integrands possess no antederivatives 
they generally possess analytic expressions for the derivative 
via the elementary calculus. As we shall see this will be 
important. In general to get the best interpolant one should 
provide good leading and trailing slopes (Q1,QN) for S(x). 

In the comparison of spline and Simpson quadrature we place 
the two techniques on equal a priori footing. This means (i) 
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that the number of data points NPTS must be odd so that 
Simpson's rule may be used; INTEG, however, only requires 

NPTS > 3, therefore let NPTS - {5,7,9,11, } for the 

comparison ; and (ii) that we compute analytically the deriva- 
t ive of the integrand so that the spline interpolant is not 
biased. Later we will discuss the problem of leading and 
trailing slopes for discrete data. In order to appraise the 
relative merits of these two techniques in the case where 
analytical methods fail, we will compare them for functions 
which do have ante derivatives. These examples provide an 
absolute reference for the quadrature determination. Since 
neither condition is violated for integrals without ante- 
derivatives the results of these comparisons should then 
follow directly for them. 

In fig. 1 are plotted the percentage absolute error for 
Simpson's 1/3 and spline quadrature for three cases. Additional 
examples are tabulated in Appendix A for a wider range of 
functions including other transcendent# is. it can be seen that 
spline quadrature in all examples has a smaller absolute error 
than does Simpson's 1/3 for a rather diverse range of functions. 
As one departs from polynomials of order - 3 for which spline 
and Simpson agree to 9 decimal degits and for 

which Simpson is exact, we see that spline quadrature is always 


PERCENT ERROR 




5 times less e ror prone for the same number of data points 
than Simpson. As the integrands develop more variability the 
absolute advantage becomes decidedly appreciable. The ^udu 
in the upper right hand corner was ’’constructed” to illustrate 
tli is fact. The integrand is shown in fig. 2. The arguments 
are driven non-linearly which gives rise to the aperiodicity 
so that fortuitous equispacing of data samples could not give 
an a_ prior advantage to Simpson. The absolute error after 61 
equispaeed data points is 2.79% vs. 12.65% for spline and 
Simpson respectively. This is ~ 10% better absolute sensitivity 
of spline quadrature versus Simpson. Though the absolute size 
of the error may vary with the interval, the 5 times more 
accurate statement is also seeln, to hold in this example. This 
greater flexibility of the quadrature is a direct manifestation 
of the global smoothness of Che interpolant in contrast to the 
jointed polynomial interpolations which are the foundation 
of the multip le quadrature rule s . 

Since spline quadrature is more accurate for a fixed tabula 
tion >f a known function whose antederivative is^ known, the 
same relative merit of spline vs Simpson quadrature is suggested 
for integrands whose antederivatives are not known, since the 
analytical existence of the antederivative in the comparison 
above was not used. 
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A PLOT OF THE INTEGRAND 
- Cos [“ exp( -|)] Sin[j exp( |-)] exp(-|-) 


Figure 2. 


A graph of the integrand udu which is asynchronously 
driven and used in the comparison in figure 1. 


• i '"liti”----™ 



GENERAL QUADRATURES 


For discrete sampling from an unknown distribution we 
seemingly have replaced the problem of numerical quadrature 
with that of determining the derivatives of the unknown distri- 
bution function so as to provide the initial and trailing slopes 
for the interpolant S(x). Numerical derivatives are a. source 
of pitfalls by themselves. The current interpolative scheme 
of Thompson is computationaly speaking, loosely connected, i.e. 
the sparse matrix involved in the linear solution for the 
spline is mostly zeroes. (Thompson, 1970 p. 5) The value of the 
leading and trailing slopes (Q1 and QN below) will propagate 
through the resulting solution for the smoothest curve which 
interpolates all the data with these values as boundary conditions. 
The loosely connected structure of the algorithum implieSj however , 
that the effects of a poor initial slope for a given set of data 
will not cause ringing more than two to three intervals away 
from the end point . To illustrate the extent to which bad 
leading and trailing slopes can affect the interpolant we now 
consider figures 3 and 4. In figure 3 we have plotted isometrically 
X = abscissa for the function y =»4 X + 5. Y in the figure is 
the percentage error of the spline interpolant vs y at each X, 
when QN = -4 and Ql = Z. Plotted here are a series of such 
error plots for Z in the range [-14,6] in two unit intervals. 
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Figure 3. Consult text page 12 



Since the spline interpolates every data point, the nodes in Y 


correspond to the x positions of the data points. The trailing 
slope in example three was obtained from the analytic form of 
y(x) and is exact. The leading slope, Ql, is varied from extreme 
to extreme to show that very poor approximations to the leading 
slope of the line propagate at most 3 data points away even 
with 10° errors in the leading slope. Figure 4 shows a 
similar result when both Ql and QN are simultaneously in error. 
The same loosely coupled result is clearly shown. 

In practice the linear approximation to the data will 
not be a disastorous mistake or be the cause of serious ringing. 
Making sure that a reasonable slope, i.e. one not inconsistent 
with the trend of the data, is a small price to pay for a 
technique otherwise free of spurious content. An example of 
such a strategy is shown in figure 5. The X and Z axes are 
the same as in figures 3 and 4, whereas the Y axis is now the 
number of equally spaced data points used in approximating the 
distribution sampled from y = sin x + 50. The leading and 
trailing slopes used for the interpolant were linear approxi- 
mations determined from the data. Note the changes in scale 
from figures 3 and 4. The errors under such a scheme are not 
very substantial. 

For discrete data a linear approximation to the data is 
a way of injecting minimal distortion into the interpolant. The 
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Figure 4. Consult text page 14 







loosely coupled nature of the interpolant implies that an 
important strategy for good interpolants over some specified 
experimental domain would be to sample outside this domain 
to allow the interpolant to settle down from the end effects 
which result from the slope approximation. 

It is thus the loosely coupled nature of the interpolant 
which allows one to say that the success of this numerical 
quadrature rule has not been undermined by the "hidden slopes'' 
problem of the underlying spline interpolant. 

FLEXIBILITY OF THE ALGORITHM: A Significant Advantage 

All of the above discussion has been from examples where 
the data is equispaced. The underlying interpolant is not 
restricted to data equispaced in abscissa. This adaptability 
of the spline interpolant makes the quadrature algorithm similarly 
flexible, for nothing in the algorithm depends on the spacing of 
the discrete data. Therefore in this algorithm is found the 
best quadrature determination regardless of data spacing consistent 
with the smooth distribution hypothesis . It should not be mis- 
construed from this last statement that this algorithm can be 
used blindly on any data set. One should make sure that the 
smooth distribution assumption is justifiable, i.e. interpolation 
is meaningful, before relying on the accuracy of this algorithm. 
Wildly fluctuating ordinates coupled with the global smoothness 
criterion for the underlying interpolant will often make the 
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quadrature result by this technique meaningless. This algorithm 

is built on an interpolation algorithm. If you would trust 

2 

the results of the smoothest C interpolant you can trust the 
results of the quadrature algorithm. 

IMPLEMENTATION 

The main rid auxiliary routines were written for use on 
the IBM 1800. For usage with more sophisticated compilers, 
however, there exist some obvious optimizations which will not 
be considered here. 

c t 'Pi- PUT T MF TMTFG ( 9 1 JM, I 9LPP , T PPMT ) 

1). The data points of the function to be integrated, 

(X, U (X) ) , should be transmitted through* COMMON as should NPTS, 
the number of data points to be interpolated and Q1 and QN, the 
initial and final slopes of the integrand. 

c nMMpM x (POO ) , U( ?0 0 ) , 9 ( POO ) , OFL ( 0 ) , n 1 * OM , UPT9 


2). If the data is tabular, decide whether linear approxi- 
mations to the solpes at X^ and are adequate. They are 

called Q1 and QN respectively. More elaborate empirical slopes 
can be devised and loaded as subroutine D(X) if desired. If 
the integrand is analytical it is desirable to load for D(X) 
the analytical function U’ (X) - D(X). 


CAL. LTMO p/\c?AMFTFP9 

T 91.. rip = 1- T NT FP WILL f pm put F A LTNFAP a PP P PX T M AT T PN 
TP O 1. AMP ON, THF T MITTAL AND FINAL 9LOPF9 PF THF IMTFP rj PLA.MT 
T 9 1 _ P p - p- TNTFF WILL LPPX FPP ANO IFF a FUNCTION SURPOUTINF 
H(V) FpP 90 MF OTHFR APP p nx T MAT T PN PP F V PL T C T T ANALYTICAL 
FPPM PF THF P)FP I \/A T I VF OF THF JN'T^CPAMO TM FVALUATION OF 
ft I A N n ON* 
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M.p. TN FTTHFP FAFF SHMF FUNCTION SIJRROUTINF D(X) MUST PR 

rnMPTLPn wtth t ntfg whcthfp dummy up pfal. it is most 

TMnnr-TANT THAT ONE SUPPLY THF PFST KNOWN VALUFS FOP 01 AMO ON. 
FOP FlIOTHFP OFT AILS CONCERNING T MF SPL I MF t OF * MASA-GRFC 
X -fto p_7o— ?F>1 * • FPL INF I MTFPPOLAT I ON ON A DIGITAL COMPUTFR • 

PV P.F.TWOWPRPM OR • SPL IMF OUADRATUOP* X-6Q2-71-200 RY 
.1 . O. SrUDDFP 


3). If matters of convergence of an Integral are 
important, or if one desires to see from which interval 
comes most of the integral's value a simple adjustment of the 
IPRNT parameter will dump the areas bounded by [X^,X^ + ^] 
with the tag i+.5 to the left. 

T PPNT = 1 I NT FO WTLL 

WRITE OUT THF PA NFL SUMS 
CONTRIBUTING TO THF INTEGRAL 

TPOMT = 2 I NTFG WILL NOT 


4). 

SI IRROl IT T NFF CALLPD 

QPI MP which FTLLS ARRAYS S AMD DFL. AND MUST BF CALLED 
PFFOPF SPLTT(J) IS USFO. 


5). 


FUNCTIONS PMPLHYFP 

SPLIT! J) WHICH PFTURNS THF I NTFPPOL.A TFD VALUE FDR 
S(X) AT THF MIOPOINT OF THF J • TH I M T F R V A L 
RFTWFFN nATA POINTS. 

D ( X ) which COMPUTFS A USFR SUPPLIED APPROXIMATION TO THE 
nFPTV/ATIVF IF RFTTFP APPROXIMATION THAN THE LINEAR ONE 
SUPPLIED IS DFSIREO. 


♦An example of such a dummy function subprogram would be 


FUNCTION D ( X ) 
D = 0 . 

RETURN 

END 


SUNK (JUT I Nf- INTFG ( SUM* I SLOP , I PR NT ) 

common x ( loo ) , u( loo ) , s ( lao ) , dfl( loo ) ♦ 01 ,om ,npts 

INTFG = INTEGRATE I c AN ALL PURPOSF NUMFR I CAL I NTFGRAT I DN 
ALGOPI Thf', THE FUNCT I (INAL PAIR? X(I),U(I) ARP COMMUNICATED 
THROUGH COMMON AS IS MPTS = NUMHFR OF DATA POINTS GRFATFR THAN 3. 
THESE VAP I AMU'S AND ARRAYS MUST RF DFFIMFD IN THE CALLING 
PROGRAM. THE APSCISf AF. * X ( I ) , Nf-irD ONLY RE DISTINCT AND MONO- 
TONIC INCREASING. N.R. THIS ALLOWS NUMERICAL INTEGRATION DP 
I R REGULARLY SPACED DP TABULATED DATA. 

THE PE TURNED VALUE SUM IS FQUAL TO THE DFFINITE INTEGRAL OF 
S ( X ) D X NVFR THF RANGF X ( 1 ) , X ( NPTS ) „ 

CALLING PARAMETERS 

I SLOP = 1- INTFG W ILL COMPUTE A LINEAR APPROXIMATION 

TO Q 1 AND ON. THE INITIAL AND FINAL. SLOPES OF THE INTFRPOLANT 
I SLOP = 2- INTFG WILL LOOK FOR AMD USE A FUNCTION SUBROUTINE 
D ( X / FOR SOME OTHER APPROXIMATION OR EXPLICIT ANALYTICAL 
FORM OF THE DERIVATIVE OF THE I NTFGPAND IN EVALUATION OF 
PH AND ON. 

N.R. IN EITHER CASE SOME FUNC T I ON . S URROUT I NE D(X) MUST BE 
COMPILED WITH I NT EG WHETHER DUMMY ‘’'OR RFAL . IT IS MOST 
IMPORTANT THAT ONE SUPPLY THF RFST KNOWN VALUES FOR Q1 AND ON. 
FOR FURTHER DETAILS CONCERNING THE SPLINE. CF. NASA-GSFC 
X -69? — 7 0-261. ‘SPLINE INTERPOLATION ON 1 A DIGITAL COMPUTER* 

MY R • F. THOMPSON OR 'SPLINE QUADRATURE* X-692-7 1-200 BY 
J.D. SCUDDFR 

IPRMT = 1 INTFG WILL 

WRITE OUT THF PA NFL SUMS 
CONTRIBUTING TO THE INTEGRAL 

IPRMT = 2 INTFG WILL NOT 
SUBROUTINES CALLED 

SPLN2*WHICH FILLS ARRAYS S AND DEL AMD MUST BE CALLFD 
BEFORE SPLIT* J)*IS USED. 

FUNCTIONS EMPLOYED 

SPLIT* J)*WHICH RETURNS THE INTERPOLATED VALUE FOR 
S(X) AT THE MIDPOINT OF THE J • TH INTFRVAL 
BETWEEN DATA POINTS. 

D ( X ) WHICH COMPUTES A USER SUPPLIED APPROXIMATION TO THE 
DERIVATIVE IF BETTER APPROXIMATION THAN THE LINEAR ONE 
SUPPLIED IS DESIRED. 

IF ( I SLOP-1 ) 13,13.14 

13 0 1 = ( U ( 2 ) — U< 1 ) ) / ( X ( 2 ) — X ( 1 ) ) 

QN = (U (NPTS) —U( NPTS— 1 ) )/( X ( NPTS > -X ( NPTS- 1 ) ) 

GO TO 15 

14 XL = X ( 1 ) 

XR = X (NPTS) 

0 1 = D ( XL ) 

ON = D(Xf$) 

15 CALL SPLN2 
SUM = 0.0 
M = NPTS- 1 

D O 1 2 I = 1 , M 
J = 1 + 1 

SUMH = DEL* J)**U(I )+4.*SPLIT< J)+U( J>)/6. 

IF(IPRNT-l) 12,9,12 
9 XARG =(I+J)/2. 

WRITE (3,100) XARG, SUMH 

100 FORMAT ( 1 2X , 'PANEL SUM OF i HE TOT A L • , / , F7 • 1 , E 1 1 .3 ) 

12 SUM = SUM + SUMH 
RETURN 
END 


: EATURES SUPPORTED 
NONPROCESS 


+ An example of such a dummy routine is seen in the text, p. 18. 
♦Copies of these routines may be found in Appendix B. 


ACKNOWLEDGEMENTS 


The author gratefully acknowledges many discussions 
with R.F. Thompson and the assistance of T.P. Carleton in 
optimizing the algorithm. 


REFERENCES 

Holladay, J.C., "A Smoothest Curve Approximation". Math Tables 
Aids Comp. 11 , (1957). 

Thompson, R.F., "Spline Interpolation on a Digital Computer", 
NASA-GSFC X -692-70-261 , 1970. 

Thompson, R.F., (Stability Analysis of Cubical Spline) to be 
published (1971) . 


APPENDIX A 


m its 



XMTE& 


Simpson $ 


CtfACT 


5 

0 


B26834E 

02 

0 

• 8 2 6 8 3 4 E 

02 

0 • 8 2 683 4 E 

0? 

02 

7 

0 


826834E 

0 2 

0 

•826834F 

02 

0 .826 8 3 4 E 

02 

9 

0 


826834E 

02 

0 

• 82 683 4 E 

02 

0. 82 683 4 E 

02 

1 1 

0 


826834E 

0 2 

0 

•82683 4 E 

02 

0 .826834F 

02 

1 3 

0 


826834E 

02 

0 

•826834E 

02 

0 .826834E 

02 

1 5 

0 


826834E 

0 2 

0 

♦ 8 2 6 8 3 4 E 

02 

0 .826834E 

02 

1 7 

0 


826834E 

02 

0 

•826834E 

0 2 

0 • 82 683 4E 

02 

1 9 

0 


826834E 

02 

0 

• 8 2 6 5 4 F 

0 2 

0 .826834E 

02 

2 1 

0 


826834E 

02 

0 

. 82 6834 E 

02 

0 .82 683 4E 

02 

23 

0 


82 6 83 4 E 

02 

0 

•826834E 

02 

0 • 8 2 6 8 3 4 E 

02 

25 

0 


826834E 

02 

0 

• 82 68 3 4 E 

02 

0 .826834E 

0 2 

27 

0 


826834E 

0 2 

0 

•826834E 

02 

0 .826834E 

02 

29 

0 


826834E 

02 

0 

• 82 6834 E 

02 

0 • 82 683 4E 

02 

31 

0 


826834E 

02 

0 

• 8 2 6 8 3 4 E 

02 

0 .8268345 

02 

33 

0 


826834E 

02 

0 

• 8 2 6 8 3 4 F 

02 

0 .826834E 

A 

V/ 

35 

0 


826834E 

02 

0 

• 826834E 

02 

0 •826834E 

0 2 

37 

0 


826834E 

02 

0 

•826834E 

02 

0.82 683 4E 

02 

5 

0 


389636E 

03 

0 

• 389636E 

03 

0 •389636E 

03 

7 

0 


389636E 

03 

0 

•389636E 

03 

0 .389636E 

03 

9 

0 


389636E 

03 

0 

• 3 8 963 6E 

03 

0 • 38 9 63 6E 

03 

1 1 

0 


389636E 

03 

0 

•389636E 

03 

0 .389636E 

03 

1 3 

0 


3 8 9 6 3 6 E 

03 

0 

•389636E 

03 

0 .389 636E 

03 

1 5 

0 


389636E 

03 

0 

•389636E 

03 

0*369 63 oE 

03 

1 7 

0 


389636E 

03 

0 

• 3 8963 6E 

03 

0 *3896365 

03 

1 9 

0 


389636E 

03 

0 

•389636E 

03 

0 • 38 963 6E 

03 

21 

0 


389636E 

03 

0 

•389636E 

03 

0 *38 96365 

0 3 

23 

0 


389636E 

03 

0 

•389636E 

03 

0 .389636E 

03 

25 

0 


389636E 

03 

0 

•389636E 

03 

0 • 389636E 

03 

27 

0 


389636E 

03 

0 

•389636E 

03 

0 .389636E 

03 

29 

0 


389636E 

03 

0 

• 3 8963 6E 

03 

0 .389636E 

03 

3 1 

0 


389636E 

03 

0 

• 38 9636E 

03 

0 .389636E 

03 

33 

0 


389636E 

03 

0 

• 3 8963 6E 

03 

0 .389636E 

03 

35 

0 


389636E 

03 

0 

•389636E 

03 

0 .389636E 

03 

37 

0 


389636E 

03 

0 

•389636E 

03 

0 .389636E 

03 

5 

0 


1 9D ( 

04 

0 

• 1 963 62 E 

04 

0 • 1 95852F 

04 

7 

0 


1 95827E 

04 

0 

• 1 95 953 E 

04 

0 • 1 95852 E 

04 

9- 

0 


195844E 

04 

0 

• 1 95884E 

04 

0 • 1 95852 E 

04 

1 1 

0 


1 95849E 

04 

0 

• 1 95865E 

04 

0 • 1 95 8 52 E 

04 

1 3 

0 


1 95851 E 

04 

0 

* 1 95858E 

04 

0 • 1 95852 E 

04 

15 

0 


1 95851 E 

0 4 

0 

. 1 95855E 

04 

0 • 1 95852E 

04 

1 7 

0 


1 95852 E 

04 

0 

• 1 9 5 8 5 4 E 

04 

0 • 1 9 5 8 5 2 E 

04 

19 

0 


1 95852 E 

04 

0 

• 1 95 853 E 

04 

0 . 1 95852 E 

04 

2 1 

0 


1 95852 E 

04 

0 

• 1 9 5 8 5 3 E 

04 

0*1 95852 E 

04 

23 

0 


1 95852E 

04 

0 

• 1 9 5 8 5 3 E 

04 

0 • 1 95852E 

04 

25 

0 


1 95852E 

04 

0 

• 1 95852 E 

04 

0 • 1 9585? E 

04 

27 

0 


1 95852E 

04 

0 

• 1 95 852 E 

04 

0 • 1 958 52 E 

04 

29 

0 


1 95852 E 

04 

0 

. 1 958 52 E 

04 

0 . 1 95852 E 

04 

31 

0 


1 95852 E 

04 

0 

• 1 95852E 

04 

0 • 1 95852E 

04 

33 

0 


1 95852 E 

04 

0 

• 1 95852 E 

04 

0*195852 E 

04 

35 

0 


1 95852 E 

04 

0 

• 1 95 852 E 

04 

0 • 1 95852 E 

04 

37 

0 


1 95 852 E 

04 

0 

• 1 95852 E 

04 

0,1 95852 E 

04 


2ir 



atr 

J * 2 


n 


i 


i 

l 

5 

t 


t 



wm 

XNTE6 


SIMPSON 

CtfACT 


5 

0 • 102347E 

05 

0 • 1 0334 9E 

05 

0 • 1 02548E 

05 

7 

0 . 1025055 

05 

0 • 1 02706E 

05 

0 . 1 02548E 

05 

9 

0 . 102535E 

05 

0 • 1 02598E 

05 

0 • 1 02548E 

05 

1 1 

0 . 102543E 

05 

0 • 1 02568E 

Ob 

0 • 1 02548E 

05 

13 

0 . 10 2 54 bE 

05 

0 • 1 02 558E 

05 

0 • 1 02548E 

05 

1 5 

0 • 102546E 

05 

0 • 1 02553F 

05 

0 • 1 02 54 HE 

05 

17 

0 . 102547E 

05 

0 . 1 0255 IE 

05 

0 • 1 02548E 

05 

19 

0 • 102547E 

05 

0 . 1 0255 0E 

05 

0 . 1 02548E 

05 

21 

0 • 102547E 

05 

0 . 1 02549E 

05 

0 • 1 02548E 

05 

23 

0 • 1 0 254 7 E 

05 

0 . 1 02549E 

05 

0 • 1 0 2 54 BE 

05 

25 

0 • 102548E 

05 

0 • 1 02548E 

0 5 

0 • 1 02548E 

05 

27 

0 • 102548E 

05 

0 • 1 02548E 

05 

0 • 1 02 54 HE 

05 

29 

0 • 102548E 

05 

0.1 02 548E 

05 

0 • 1 02548E 

05 

31 

0 • 102548E 

05 

0 . 1 02 54 BE 

05 

0 • 1 02548E 

05 

33 

0 • 10 2.548E 

0 5 

0 . 1 02548E 

05 

0 • 1 02548E 

05 

35 

0 • 102548E 

05 

0 . 1 02548E 

05 

0 • 1 0 2 54 BE 

05 

37 

0 . 102548E 

05 

0 • 1 02548E 

05 

0 . 1 02548E 

05 

5 

0 • 54978 f\z 

05 

0 .561 900E 

05 

0 . 5522 82 E 

05 

7 

0 .551 786E 

05 

0 .55423 IE 

05 

0 .552282L 

05 

9 

0 ,5521 25E 

05 

0 .552904E 

05 

0 . 552282 E 

05 

1 1 

0.55221 7E 

05 

0 .552538E 

05 

0 • 552.2 82 E 

05 

1 3 

0 • 552251 E 

05 

0 .552405E 

05 

0 • 552 2 82 E 

05 

15 

0 .552265E 

05 

0 .552349E 

05 

0 .552282E 

05 

17 

0 • 552272E 

05 

0.55232 IE 

05 

0 • 552282 E 

05 

19 

0 .552276E 

05 

0 .552306F 

05 

0 . 552 2 82 E 

05 

21 

0 .552278E 

05 

0 .552298E 

05 

0 .5522.82E 

05 

23 

0 .552279E 

05 

0 .552293E 

05 

0 .5522 82E 

05 

25 

0 .552280E 

05 

0 .552289E 

05 

0 • 552 2 82 E 

05 

27 

0 .552280E 

05 

0 • 5522 87E 

05 

0 .552282E 

05 

29 

0 .552281 E 

05 

0 .552286E 

05 

0 • 5522 82 E 

05 

31 

0 .552281 E 

05 

0 .552285F. 

05 

0 .552382E 

05 

33 

0 . 55228 1 E 

05 

0 .552284E 

05 

0 . 5522 82 E 

05 

35 

0 .552281 E 

05 

0 .552284E 

05 

0 . 5522 82 E 

05 

37 

0 . 552281 E 

05 

0 .5522 83E 

05 

0 .552282E 

05 

DRp 







5 

0 .30091 4E 

06 

0 .3137 14E 

06 

U • 8 03 68^ r: 

06 

7 

0 . 303090E 

06 

0 .305732E 

06 

0 .303632E 

06 

9 

0 . 30 346 0E 

06 

0 .304309E 

06 

0 .303632E 

06 

1 1 

0 . 30 3 562 E 

06 

0 .3039 12E 

06 

0.30 3632 E 

06 

13 

0 .303598E 

06 

0 .303768E 

06 

0 • 3 03 632 E 

06 

1 5 

0 .30361 4E 

06 

0 .303706F 

06 

0.30 3632 E 

06 

17 

0 .303622E 

06 

0 .303675E 

06 

0 .303632E 

06 

19 

0 .303626E 

06 

0 .303659E 

06 

0 .3 0 3 632 E 

06 

21 

0 .303628E 

06 

0 .303650E 

06 

0 .303632E 

06 

23 

0.303629E 

06 

0 . 303645E 

06 

0 .303632E 

06 

25 

0 .303630E 

06 

0 .30364 IE 

06 

0 .303632E 

06 

27 

0 . 303631 E 

06 

0 .303639E 

06 

0 . 3 03 632 E 

06 

29 

0 .303631 E 

06 

0 .3 03 63 7E 

06 

0 .303632E 

06 

31 

0 .303632E 

06 

0 .303636E 

06 

0 .303632E 

06 

33 

G • 30 3632E 

06 

0 .303635E 

06 

0 .303632E 

06 

35 

0 . 30 3632 E 

06 

0 .303635E 

06 

0 .303632E 

06 

37 

0 .303632E 

06 

0 . 3 03 634 E 

06 

0 .3 03 632 E 

06 

5 

0 • 1 66 88 0E 

07 

0.1791 ' 

07 

0 . 1 6958 OE 

07 

7 

0 • 1 69038E 

07 

0.171 fo3c£ 

07 

0 • 1 69580E 

07 

9 

0 • 169408E 

07 

0 . 1 70250E 

07 

0 • 1 6958 OE 

07 

1 1 

0 . 169509E 

07 

0 • 1 69858E 

07 

0 • 1 69580E 

07 

1 3 

0 . 169546E 

07 

0 . 1 697 15E 

07 

0 . 1 6958 OE 

07 

15 

0 • 1 69562 E 

07 

0.1 69653E 

07 

0 . 1 69580E 

07 

17 

0 • 169569E 

07 

0. 169623E 

07 

0 • 1 6958 0 E 

07 

19 

0 • 1 69573E 

07 

0 • 1 69607E 

07 

0 • 1 695 8 0 E 

07 

21 

0 • 169576E 

07 

0 • 1 69598E 

07 

0 • 1 6958 OE 

07 

23 

0 . 169577E 

07 

0 . 1 69592E 

07 

0 • 1 69580E 

07 

25 

0. 169578E 

07 

0 . 1 69589E 

07 

0 • 1 69580E 

07 

27 

0 • 169579E 

07 

0 • 1 69586E 

07 

0 . 1 69580E 

07 

29 

0. 169579E 

07 

0 . 169585E 

07 

0 • 1 6958 OE 

07 

31 

0 • 1 69579E 

07 

0 • 1 69584E 

07 

0.1 69580E 

0 7 

33 

0. 169579E 

07 

0 • 1 69583E 

07 

O'. 1 6958 OE 

07 

35 

0 • 1 69580E 

07 

0 • 1 69582E 

07 

0 • 1 69580E 

07 

37 

0. 169580E 

07 

0 • 1 69582 E 

07 

0 • 1 69580E 

07 


23 * 




- ' r • ' - x ~ (iniiiii V i - ~ - 


2tr 

jit s d< 


o 




0 






Nfts 

TNTE6 


SintpsoW 

% 

E*ACT 


5 

0 .9335535 

0 7 

0 . 1 042275 

0 8 

0 .9389565 

0 7 

7 

0 .9538755 

07 

0.9776905 

0 7 

0 • 9589365 

07 

9 

0 .9573355' 

07 

0.9631325 

0 7 

0.9589365 

07 

1 1 

0 .9582895 

07 

0.961 3465 

0 7 

0 • 9 3 895 6F 

0 7 

1 3 

0 .9585348. 

0 7 

0 .96 02 19F 

07 

0 . 9 3 8 9 3 6 F 

0 7 

1 5 

0 .938781 5 

07 

0 . Q 3 9 6 4 2 F 

07 

0 . 9 3 8 93 65 

0 7 

1 7 

0 .9 388335 

0 7 

0 .9 393605 

07 

O' . 9 3 8 9 3 65 

0 7 

1 9 

0 .9 588 92 5 

0 7 

0.95 92 0 8 E 

07 

0 .9389365 

07 

2 1 

0.9389145 

0 7 

0 .939 1 22 E 

0 7 

0 .93893 65 

0 7 

23 

0 .9 38 92 75 

07 

0 • 9 5 9 o 6 9 E 

07 

0 . 9389365 

07 

25 

0 .9389355 

0 7 

0 .959 03 6E 

07 

0 . 9389365 

0 7 

27 

0.938941 5 

0 7 

0 .9390 1 4 E 

0 7 

0 .93893 65 

07 

29 

0 .9389455 

0 7 

0 .9589995 

07 

0.9369365 

0 7 

3 1 

0 .9 58 94 75 

07 

0 .958989E 

07 

0.93695 65 

07 

33 

0 .9589495 

0 7 

0 .95898 IE 

07 

0 • 9389365 

07 

35 

0.9589515 

07 

0 . 958976E 

07 

0 .93 8 93 65 

07 

37 

k ^~EF 

0 .9589525 

07 

0 .95897 IE 

07 

0 .9389365 

07 

1 

5 

0.52 3643E 

08 

0 .61 61 92E 

08 

0 . 5477345 

08 

7 

0 .543232E 

08 

0 • 5 63 9 1 6E 

08 

0 .3477 345 

OH 

9 

0 • 54630 6E 

08 

0.5531 93 E 

08 

0 . 5477345 

08 

1 1 

0 .5471 57E 

08 

0 .55 0 0 4 6E 

08 

0 .5477345 

08 

1 3 

0 . 5474655 

08 

0 .548876E 

08 

0.5477545 

08 

1 5 

0 • 547598E 

08 

0 .5483 6 5E 

08 

0 .5477545 

08 

1 7 

0 . 54 7662 E 

08 

0.5481 15E 

08 

0 . 54 775 45 

08 

19 

0 • 547697E 

0 8 

0 .54798 OE 

08 

0 .54 77 345 

08 

2 1 

0 . 5477165 

08 

0 .547903E 

08 

0 • 5477345 

08 

23 

0 • 547 72 BE 

08 

0 .547856F 

0 8 

0 .5477 345 

08 

25 

0 .547736E 

08 

C .5478266 

08 

0 .547754E 

08 

27 

0 .54 774 IE 

08 

0 .547806 E 

08 

0 .547754E 

08 

29 

0 • 547744E 

08 

0.5477935 

08 

0.5477345 

08 

3 1 

0 . 54 7 74 7E 

08 

0 .5477835 

08 

0 . 5477545 

08 

33 

0 . 547748E 

08 

0 .54 77 775 

08 

0 . 547734E 

08 

35 

0 .547750E 

0 8 

0 .547772E 

08 

0 .547734E 

0 8 

37 

0 . 3477305 

08 

0 .547768E 

08 

0 . 54 77 54 5 

OB 




5 

0 .7 1 52 31 669 E 

00 

7 

f) • 7 1 5233 658 E 

00 

9 

0 .7 1 52 3 3993 F 

00 

1 1 

0.71 5234084E 

00 

1 3 

0.71 52341 1 7F 

00 

15 

0 .7 1 52341 30 F 

00 

17 

0 .7 1 52341 36F 

00 

1 9 

0 .71 52341 39 F 

00 

21 

0 .71 52341 42 F 

00 

23 

0.71 5234144F 

00 

25 

0.71 5234143F 

00 

27 

0.71 52 34 1 44F 

00 

29 

0.71 52341 43E 

00 

31 

0.71 5234143F 

00 

33 

0 .7 1 52341 44E 

00 

35 

0.71 52341 43E 

00 

37 

0.71 52341 44E 

00 


0.715244137E 00 
0.71 52361 UE 00 
0.715234767E 00 
0 • 7 1 52 3 4 4 0 2 F 00 
0.715234270E 00 
0*71 52342 1 3E 00 
0*7152341 86E 00 
0*71 52341 7 0E 00 
0* 7152341 62 E 00 
0*71 52341 57E 00 
0*71 5234 1 54E 00 
0*71 52 341 5 0 E 00 
0*71 5234 1 49E 00 
0*71 52341 49E 00 
0.7152341455 00 
0*71 52 3414 7E 00 
0.715234146E 00 


FXACT 


0 .71 52-741 45F 00 


1 


7/4 

xdx 


24 


i n ,■ MMMliMiir^iii«ntMMTiii>iiil ia i a^ r m - iiiiiiHiniaflMr- • i r - ' jrr1ri' , -i - 'lir' 


NPTS SPLINE 


Simpson 


5 

0 


29289261 3F 

7 

0 


292893099F 

9 

0 


29 289 3 1 8 OF 

1 1 

0 


29 289 320 3F 

1 3 

0 


2928932 10E 

15 

0 


29289321 4F 

17 

0 


29289321 5E 

19 

0 


29289321 6E 

21 

0 


2928932 1 7E 

23 

0 


29289321 7F 

25 

0 


2928932 1 7E 

27 

0 


29289321 7F 

29 

0 


29289321 7F 

31 

0 


29289321 7F 

33 

0 


2928932 1 7F 

35 

0 


2 9289321 7E 

37 

0 


29289321 7E 


00 

0 


292895648F 

00 

0 


292 893697E 

00 

0 


292R93369E 

00 

0 


2 92 893280E 

00 

0 


292893248F 

00 

0 


292 8 9323 4E 

00 

0 


29289322 7E 

00 

0 


292893224E 

00 

0 


292893222E 

00 

0 


292 89322 CE 

00 

0 


29289322 OF 

00 

0 


29289321 9E 

00 

0 


29289321 8F 

00 

0 


29289321 9E 

00 

0 


2928932 1 8F 

00 

0 


29289321 8E 

00 

0 


2928932 1 8F 


00 
00 
00 
00 
00 
00 
00 
0 0 

00 0.2928932 1 9E 00 



5 

0 


1 57 1 60 68 4 E -07 

0 


1 46291 R07E-07 

7 

0 


21 536R345F-07 

0 


21 374 85 8 6E — 0 7 

9 

— 0 


183354131 F — 0 8 

-0 


79241 3960E-08 

1 1 

0 


366708 2 63F — 0 8 

0 


1 1 21 57052E— 08 

13 

0 


19615981 7F-07 

0 


1 85 70 9322 E— 07 

15 

0 


4249 1 5925F-08 

0 


247302818F— 08 

17 

-0 


416 184775E-08 

— 0 


512021 327E-08 

19 

-0 


17 17 12599E-08 

0 


647476705E-08 

21 

0 


28885551 7F-08 

“0 


853368880E— 09 

23 

0 


8207280 19F-08 

-0 


74 6974837E- 08 

25 

0 


226827978 F— 08 

-0 


692 8542 56E— 08 

27 

-0 


982936399E— 09 

0 


3 9 57380 96E — 08 

29 

0 


13551 47 1 0 E—09 

0 


44061 6992E-08 

31 

0 


41 1 09 1 605F— 09 

0 


3072 1 2796E-08 

33 

-0 


6 00 26650 3 E—09 

0 


472 40 062 9E — 09 

35 

-0 


26 19344 74E— 08 

-0 


573693364E— 09 

37 

0 


97861 6299F-09 

0 


32780201 4E— 08 


0.0 


J CoC(l00*>d* 


5 

0 


2 1 444351 3fc 

DO 

7 

0 


2 1 45698 39E 

00 

9 

0 


214 59 1 62 7F 

00 

1 1 

0 


21 4597637F 

00 

1 3 

0 


21 4599806F 

00 

15 

0 


21 4600739F 

00 

17 

0 


21460 1 192E 

00 

19 

0 


21460 1 434F 

00 

21 

0 


21460 1572E 

00 

23 

0 


21460 1 655F 

00 

25 

0 


21 460 1708F 

00 

27 

0 


21460 1 743E 

00 

29 

0 


21 460 1766E 

00 

31 

0 


2 1 460 1 7 8 3F 

00 

33 

0 


21 460 1795F 

00 

35 

0 


2 1 460 1 8 04E 

00 

37 

0 


2 1 4 60 1 8 1 0 E 

00 


O • 2 1 bl bO MOt 00 
0.214721053E 00 
0.214640998E 00 
0.21461R174E 00 
0.214609796E 00 
0.214606159E 00 
0*214 604381 E 00 
0*21 46 03 42 9E 00 
0.214602883E 00 
0 *21 4602552E 00 
0.214602342E 00 
0*21 4602204E 00 
0*2146021 09E 00 
0 *21 4602043E 00 
0.214601996E 00 
0*21 4601 962E 00 
0.214601 93 6E 00 


0.214601834E 00 



5 

0 i 

. 1 19327758E 

0 1 

7 

0 . 

. 1 1 9 327956E 

01 

9 

0« 

► 1 19327989E 

0 1 

1 1 

0 . 

, 1 19327998E 

0 1 

13 

0 , 

. 1 1932800 IF 

0 1 

15 

0 . 

. 1 19 328003E 

01 

17 

0 « 

» 1 19328003F 

0 1 

19 

0 . 

. 1 19328004F 

01 

21 

0 « 

► 1 19328004F 

0 1 

23 

0 . 

, 1 19328004F 

0 1 

25 

0 i 

► 1 19328004E 

0 1 

27 

0 . 

. 1 19328004F 

01 

29 

0 , 

. 1 19328004E 

0 1 

31 

0 . 

. 1 19328004E 

0 1 

33 

0 , 

. 1 19328004F 

0 1 

35 

0 . 

, 1 1932R004F 

0 1 

37 

0 i 

► 1 19328004F 

0 1 


0 

0 

0 

0 

O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


119328985E 01 
119328199E 01 
119328066E 01 
1 1 932 8 03 OE 01 
1 1 93280 1 7E 01 
1 1 93280 1 1 E 01 
1 1 9328008E 01 
1 1 9328007E 01 
1 1 93280 06E 01 
1 1 9328005E 01 
1 1 9328005E 01 
1 1 9328005E 01 
1 1 9328005E 01 
1 1 932 8004E 01 
1 1 °3280 04E 01 
119328004E 01 
1 1 93280 04E 01 


0*1 19328004E 01 



um 

SPLINE 



SlMPSOlJ *(s 



€XACT 


6 

U . 1 26487651 E 

04 

0 

• 1 352 82 39 1 E 

04 

0 

• 1 286966 09E 

04 

7 

0 • 1 26228003E 

04 

0 

• 126382290E 

04 

0 

. 1 26 25985 7E 

04 

9 

0 • 1 26249733E 

04 

0 

• 1 262 99443E 

04 

0 

. 1 26259857E 

04 

1 1 

0 . 1 262 5570 IE 

04 

0 

• 1 26276236E 

04 

0 

. 1 2625985 7E 

04 

1 3 

0 • 1 26 257 8 50 E 

04 

0 

• 126267800E 

04 

0 

• 1 26259857E 

04 

1 5 

0 • 126258773E 

04 

0 

• 1 262641 58E 

04 

0 

. 1 2625985 7E 

04 

1 7 

0 • 1 2625922 1 E 

04 

0 

• 1 26262384E 

04 

0 

• 1 26259857E 

04 


0 . 1 262 59 4 59E 

04 

0 

• 126261 437E 

04 

0 

. 1 2625985 7E 

04 

2 1 

0 • 1 26 2 59 59 6F 

04 

0 

• 1 2626 0894 E 

04 

0 

• 1 26259857E 

04 

23 

0 . 126259678E 

04 

0 

• 1 26260566E 

04 

0 

. 1 2625985 7E 

04 

25 

0 . 1 26 2 597 30 E 

04 

0 

• 1 26260357E 

04 

0 

• 1 26259857E 

04 

27 

0 * 126259765E 

04 

0 

• 126260220E 

04 

0 

. 1 2 6 25985 7E 

04 

29 

0 . 1 26259788E 

04 

0 

• 126260127E 

0 4 

0 

. 1 26 259857E 

04 

3 1 

0 • 1 262 59 8 0 5E 

04 

0 

. 1 26260062E 

04 

0 

• 1 2625985 7E 

04 

3 3 

0. 1 262598 1 6E 

04 

0 

«, 1 2626001 5E 

04 

0 

. 126259857E 

04 

35 

0 . 1 26259825E 

04 

0 

• 1 2 62 5998 1 E 

04 

0 

• 1 2 6 25985 7E 

04 

37 

0. 126 2 598 3 1 E 

04 

0 

• 1 262 5 9955F 

04 

0 

• 126259857E 

04 

5 

0 • ltcoiU2531E 

04 

0 

. 1 26845534E 

04 

0 

• jl 8626 1 7 88 1 

04 

7 

0 . 1 26229933E 

04 

0 

• 1 2 63 842 2 3E 

04 

0 

• 1 26261788E 

04 

9 

0 . 1 26251 664 E 

04 

0 

. 12630 1 375E 

04 

0 

. 12626 1788E 

04 

1 1 

0 • 1 262 57 632 E 

04 

0 

• 1262781 68E 

04 

0 

. 1 26261788E 

04 

1 3 

0 . 126259782E 

04 

0 

• 126269731 E 

04 

0 

• 1 26 26 1 7 88E 

04 

1 5 

0 • 1 26260704E 

04 

0 

• 12626O089E 

04 

0 

• 1 26261788E 

04 

1 7 

0 . 126261 1 52 E 

04 

0 

• 12626431 5E 

04 

0 

. 12626 1788E 

04 

19 

0 . 126261 391E 

04 

0 

. 126263368E 

04 

0 

• 1 2626 1788E 

04 

2 1 

0 . 1 26261 52 7E 

04 

0 

. 1 262 62 82 5E 

04 

0 

• 1 2626 1788E 

04 

23 

0 . 1 26261 6 1 0 E 

04 

0 

• 1 262 62497E 

04 

0 

. 1 26261788E 

0 4 

2 5 

0.1262 61 662 F. 

04 

0 

• 1 26262289E 

04 

0 

• 1 2626 1788b 

04 

27 

0 . 126261 697E 

04 

0 

. 1262621 52E 

04 

0 

• 1 2626 1788E 

04 

29 

0 . 1 26261720E 

04 

0 

. 1 26262 058E 

04 

0 

• 12626 1788E 

04 

31 

0 . 1 26261 736E 

04 

0 

• 126261 993E 

04 

0 

. 1 2626 1 7 88E 

04 

33 

0 . 1 262 617 48 E 

04 

0 

• 126261 946E 

0 4 

0 

• 1 2626 1788E 

04 

35 

0 . 1262617 56E 

04 

0 

• 126261 91 2E 

04 

0 

• 1 26261788E 

04 

37 

0 • 1 262 617 63 E 

04 

0 

• 1 26261 887E 

04 

0 

• 1 2626 1 788E 

04 

5 

0 .3926610 1 5E 

0 1 

0 

.392 6581 2 6E 

0 1 

0 

. 39266 0 28 3E 

0 1 

*7 

f 

0 • 39 2660 4 3 5 E 

0 1 

0 

. 3926597 6 1 E 

01 

0 

•392660283E 

0 1 

9 

0 .392660332E 

0 1 

0 

.392 660 1 04E 

0 1 

0 

•392660283E 

0 1 

1 1 

0 .392660 303E 

0 1 

0 

. 3926602 07E 

01 

0 

•39 266028 3E 

0 1 

1 3 

0 . 392 660293E 

0 1 

0 

. 3 92 660245E 

0 1 

0 

. 39 266 0 28 3E 

0 1 

1 5 

0 . 392660 288E 

0 1 

0 

• 3926602 62E 

01 

0 

• 39266028 3E 

0 1 

1 7 

0 .392660285E 

0 1 

0 

. 392 6602 7 1 E 

0 1 

0 

. 39266 0 28 3E 

0 1 

1 9 

0 . 392660 284E 

0 1 

0 

• 392660275E 

01 

0 

• 39266028 3E 

0 1 

2 1 

0 .392660283E 

01 

0 

. 392 66 0 2 77E 

0 1 

0 

• 392660283E 

0 1 

23 

0 • 392660 283E 

0 1 

0 

. 3926602 7 9E 

01 

0 

•392660283E 

0 1 

25 

0 .392660282E 

0 1 

0 

. 392 660279E 

0 1 

0 

. 39 266 0 28 3E 

01 

27 

0 • 392660 282E 

0 1 

0 

•392660280E 

01 

0 

• 39 266028 3E 

0 1 

29 

0 .392660283E 

0 1 

0 

. 392 660280E 

0 1 

0 

. 3 9 266 0 28 3E 

01 

31 

0.392660282E 

0 1 

0 

•392660280E 

01 

0 

• 3 9 266028 3E 

0 1 

33 

0 .392660282E 

0 1 

0 

.392660281 E 

0 1 

0 

• 392660283E 

0 1 

35 

0 • 3926 60 28 1 E 

0 1 

0 

• 39266028 1 E 

01 

0 

• 3 9 2660 28 3E 

0 1 

37 

0 • 39266028 1 E 

0 1 

0 

.392660281 E 

0 1 

0 

. 392660283E 

0 1 

5 

0 .3927371 57E 

0 1 

0 

• 392 74 0057E 

0 1 

0 

. 3 927 378 9 4E 

01 

7 

0 • 3927 377 41 E 

0 1 

0 

.39273841 8E 

01 

0 

• 3 9 27 3789 4E 

0 1 

9 

0 .392737844E 

0 1 

0 

. 392 738073E 

01 

0 

. 3927 37894E 

01 

1 1 

0 .392737873E 

0 1 

0 

. 392737969E 

01 

0 

• 39 27 3789 4E 

0 1 

1 3 

0 .392737883E 

0 1 

0 

.392737931 E 

0 1 

0 

. 3927 378 94E 

01 

15 

0 .392737887E 

0 1 

0 

.39273791 3E 

01 

0 

• 3 9 27 37894E 

0 1 

17 

0 .392737889E 

0 1 

0 

. 392 737905E 

01 

0 

• 39 27 37 8 9 4E 

0 1 

19 

0 . 39 27 37 8 90 E 

0 1 

0 

. 392737900E 

01 

0 

• 39 27 3789 4E 

0 1 

2 1 

0 .392737891 E 

0 1 

0 

. 392 73 7898E 

0 1 

0 

.39273789 4E 

01 

23 

0 .392737891E 

0 1 

0 

• 392 7378 96E 

0 1 

0 

• 3927 37894E 

0 1 

25 

0 .392737891 E 

0 1 

0 

. 392 737895E 

01 

0 

•392737894E 

01 

27 

0 .392737892E 

0 1 

0 

. 392737894E 

01 

0 

• 3 9 27 3789 4E 

o ; 

29 

0.392737891E 

0 1 

0 

. 392 73 7893E 

0 1 

0 

• 3927 37894E 

01 

31 

0 .392737891E 

0 1 

0 

•392737893E 

01 

0 

• 3 9 27 37894E 

0 1 

33 

0 .392737891 E 

0 1 

0 

. 392 73 7893E 

01 

0 

• 3 9 27 37 8 9 4E 

01 

35 

0 .392737891E 

o 1 

0 

. 392737 892E 

01 

0 

• 3 9 27 37 89 4E 

0 1 

37 

0 .392737892E 

0 1 

0 

. 392 73 7 8 92 E 

01 

0 

• 3927 3789 4E 

01 
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APPENDIX B 


SUBROUTINE SPLN2 
DIMENSION A(100),V(100) 

COMMON X< 100), U(100),S( 100), DEL(100),Q1,QN, NPTS 
UNEQUALLY SPACED DATA 

THIS PROGRAM COMPUTES THE SECOND DERIVATIVE, S(X), OF THE 
CUBIC SPLINE, SPLIN(X), WHICH INTERPOLATES THE NPTS OF 
ARBITRARILY SPACED DATA <X,U). 

THE COMMON STATEMENT PROVIDES COMMUNICATION WITH THE 
FOLLOWING FUNCTION SUBROUTINES 

FUNCTION 

SUBROUT I NE DESCRI PTI ON 

SPLIN(X) FOR THE VALUE X, SPLIN RETURNS THE VALUE OF THE 
CUBIC SPLINE WHICH INTERPOLATES THE NPTS 
DATA POINTS <X,U> 

DSPLN(X) RETURNS THE DERIVATIVE OF THE SPLINE AT X. 


N-NPTS 

I F(N-3) 5,5,1 

1 I F(N-100) 6,6,5 

COMPUTE DEL AND V 

6 DEL(2)-X(2)-X<1) 

V(l)-6.0*(( (U(2)-U(l) )/DEL(2 ) )-Ql ) 

Nl-N-1 

DO 2 I -2,N1 

DEL(l+l) a X(l+l)-X(l) 

2 V(l)»((U(l-l)/DEL(l))-U(l)*((1.0/DEL(l))+(1.0/DEL(|4l))) 
C +<U(U1)/DEL(I^1)))*6.0 

V(N)«(QN+CU(N1)-U(N))/DEL(N) )*6.0 

GAUSSIAN ELIMINATION AND AUGMENTATION 

A(1)-2.0*DEL(2) 

A( 2 )-l . 5*DEL(2 )+2 .0*DEL( 3) 

V(2)-V<2)-0.5*V<1) 

DO 3 1-3, N1 
C-DELC I )/A( 1-1) 

A(I)-2.0*(DEL(I )*DEL( 1+1) )-C*DEL( I ) 

VC I )-V( I )-C*VC 1-1) 

3 CONTINUE 
C-DEL(N)/A(N1 ) 

A(N)-2.0*DEL(N)-C*DEL(N) 

V(N)-V(N)-C*V<N1) 

BACK SUBSTITUTION 
S(N)-V(N)/A(N) 

DO 4 J-1,N1 
l-N-J 

4 S(I)-(V(I)-DEL(I+1)*S(I+1))/A(I) 

5 RETURN 
END 
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APPENDIX B 


SUBROUTINE SPLN2 
DIMENSION A(X00) / V(100) 

COMMON X(100),U(100),S(100),DEL(100 ) ,Q1,QN, NPTS 
UNEQUALLY SPACED DATA 

THIS PROGRAM COMPUTES THE SECOND DERIVATIVE, S(X), OF THE 
CUBIC SPLINE, SPLIN(X), WHICH INTERPOLATES THE NPTS OF 
ARBITRARILY SPACED DATA (X,U). 

THE COMMON STATEMENT PROVIDES COMMUNICATION WITH THE 
FOLLOWING FUNCTION SUBkOUTINES 

FUNCTION 

SUBROUTINE DESCRIPTION 

SPLIN(X) FOR THE VALUE X, SPLIN RETURNS THE VALUE OF THE 
CUBIC SPLINE WHICH INTERPOLATES THE NPTS 
DATA POINTS <X,U> 

DSPLN(X) RETURNS THE DERIVATIVE OF THE SPLINE AT X. 


N-NPTS 

I F (N-3 ) 5,5,1 

1 IF(N-IOO) 6,6,5 

COMPUTE DEL AND V 

6 DEL( 2 ) "X ( 2 ) - X ( 1 ) 

V(1)-6.0*(((U(2)-U(1))/DEL(2))-Q1 ) 

Nl-N-1 

DO 2 1-2, N1 

DEL( l*l)-X( l+D-XC I ) 

2 V( I )-((U( l-l)/DEL( I ))-U( I )*((1.0/DEL( I ))+(1.0/DEL( !♦!))) 
C ♦(U(I+1)/DEL(I^1)))*6.0 

V(N)-(QN^(U(N1)-U(N))/DEL(N) )*6.0 

GAUSSIAN ELIMINATION AND AUGMENTATION 

A( 1 )-2 • 0*DEL( 2 ) 

A(2)-1.5*DEL(2)**2.0*DEL<3) 

V(2)-V(2)-0.5*V(1) 

DO 3 1-3, N1 
C-DELC I )/A( I -1 ) 

A( I )-2 .0*(DEL( I )^DEL( I ^1) )-C*DEL( I ) 

V(I)-V(I)-C*V(I-1) 

3 CONTINUE 
C-DEL(N)/A(N1 ) 

A(N)-2.0*DEL(N)-C*DEL(N) 

V(N)"V(N)-C*V(N1) 

BACK SUBSTITUTION 
S(N)-V(N)/A(N) 

DO 4 J-l, N1 
l-N-J 

4 S(I)-(V(I)-DEL(I«1)*S(I«1))/A(I) 

5 RETURN 
END 
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FUNCTION SPLIT(N) 


THIS SUBROUTINE IS MEANT TO BE USED IN CONJUNCTION WITH 
SPLN2 AND I NT FG AND RETURNS THE VALUE OF THE CUBICAL 
SPLINE INTERPOLANT AT THE MIDPOINT OF THE N'TH INTERVAL. 
BETWEEN DATA POINTS. 

COMMON X ( 100 ) *U< 100 ) . S( 100 > . DEL ( 100 ) *01 ,QN.NPTS 
D = DEL ( N ) / 2 . 0 

SPLIT = (U (N) +U( N-l ) ) /2.0-D*D* ( S< N) +S< N- 1 ) )/4.0 
RETURN 
END 


FEATURES SUPPORTED 


